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Non-equilibrium spin transport in a magnetic tunnel junction comprising a single magnetic layer 
in the presence of interfacial spin-orbit interaction (SOI) is studied theoretically. The interfacial 
SOI generates a spin torque of the form T = r||M x (z x M) -I- T±z x M, even in the absence of an 
external spin polarizer. For thick and large tunnel barriers, the torque reduces to the perpendicular 
component, T±, which can be electrically tuned by applying a voltage across the insulator. In the 
limit of thin and low tunnel barriers, the in- plane torque T\\ emerges, proportional to the tunneling 
current density. Experimental implications on magnetic devices are discussed. 



INTRODUCTION 

Current-driven magnetization dynamics (switching 
and excitations) using Spin Transfer Torque [li (STT) 
has essentially been observed in inhomogeneous magnetic 
structures such as spin-valves, magnetic tunnel junctions 
and magnetic domain walls Indeed, in such struc- 
tures an external polarizer or, equivalently, an inhomo- 
geneous spin texture is needed to allow for angular mo- 
mentum transfer from the spin momentum of itinerant 
carriers to the localized magnetic moments. Experimen- 
tally, this implies either complex magnetization dynamics 
(complex domain wall motion and transformations Q ) or 
advanced design involving multilayered structures such 
as synthetic antiferromagnets Q to stabilize the exter- 
nal polarizer. Therefore, controlling the magnetization 
direction of a single ferromagnet without the need of an 
external polarizer or magnetic texture would provide sig- 
nificant advantage for novel spin devices. 

Among the different methods available, magnetoelec- 
tric effect [1], multiferroic and magnetostrictive Q 
materials have been proposed. In Ga(Mn,As), electrical 
control of magnetic anisotropy through the voltage de- 
pendence of hole density has been demonstrated [9| and 
in rare-earth/transition metals compounds laser-induced 
magnetization reversal through inverse Faraday-effect 
has been achieved [l^. Besides all these mechanisms, 
the electrical control of interfacial magnetic anisotropy 
offers an elegant and promising solution 
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Indeed, recent experiments have demonstrated the pos- 
sibility to manipulate the interfacial magnetic anisotropy 
of a thin ferromagnet by apply ing a bias voltage through 
an adjacent insulator [ll|, [l^- Voltage-driven Magneto 
Crystalline Anisotropy (MCA) change has been observed 
for large bias voltages [l^ . This observation has been ex- 
plained in terms of voltage-controlled filling of electron 
orbitals at the interface between the insulator and the 
ferromagnet In this configuration, the bias volt- 

age generates a non-equilibrium Perpendicular Magnetic 
Anisotropy (PMA) that can in principle switch the mag- 
netization direction from in-plane to out-of-plane and 
vice-versa. Notice that this voltage-driven PMA acts as 
an effective field that can not excite self-sustained mag- 



netization precessions. 

Alternatively, a recent prediction suggested that 
the spin-orbit interaction (SOI) present in a asymmet- 
rically designed thin single ferromagnetic layer could be 
used to generate current-driven effective magnetic field 
and to manipulate the local magnetization. This predic- 
tion has been experimentally confirmed recently |15l. Il6| 
and followed by a number of theoretical investigations 
[l7| . Similarly to the voltage-driven PMA mentioned 
above, one of the drawbacks of the SOI-induced torque 



(SOI-ST) derived in Ref. 14 1 is that the torque acts as 
an effective magnetic field 15,[l6| that does not compete 



with the Gilbert damping (contrary to STT) so that no 

self-sustained precessions can be excited. 
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FIG. 1. (Color online) Schematics of a Semi-Magnetic Tunnel 
Junction comprising a single ferromagnet and a Rashba-type 
SOI at the interface between the ferromagnet and the insula- 
tor. 

In the present paper, we study the non-equilibrium 
spin transport in a Semi-Magnetic Tunnel Junction 
(SMTJ), comprising a single ferromagnetic layer as well 
as interfacial SOI at the interface between the ferromag- 
net and the tunnel barrier. At equilibrium (zero bias), 
the interfacial spin-orbit coupling is responsible for a 
large PMA. In the case of thick and large barrier, when 
tunneling conductance vanishes, the amplitude of PMA 
can be tuned by a gate voltage, yielding voltage- controlled 
PMA, as experimentally demonstrated in Ref. 
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the non-equilibrium regime and thin barrier limit, the in- 
terfacial SOI generates an additional current- driven spin 
torque on the ferromagnetic layer that can induce either 
magnetization switching or self-sustained magnetic pre- 
cessions lisl. 
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FIG. 2. (Color online) Angular dependence of the TAMR 
and perpendicular torque at zero bias for q_r = 5eV.A^ and 
d = Inm. The other parameters are given in the text. 



THEORETICAL MODEL 

SMTJs are well known to display tunneling anisotropic 
magnetoresistance 1^21| (TAMR) : the resistance of the 
junction depends on the orientation of the magnetization 
compared to the normal of the layers (angle 9 in Fig. 
1). This effect has been observed in a wide range of 
F/I interfaces, where F is a ferromagnet (Fe, Ni, Co, 
GaMnAs etc.) and I is an insulator (AlOx, MgO, GaAs, 
etc.) [3" 21 1 . Interestingly, TAMR is interpreted as an 
effect of interfacial spin-orbit coupling, such as Rashba 
and Dresselhaus SOI fHlli: 



ayky)]S{z) 



Hsoi = [aR{axky - aykx)S{z) + aoicTxkx 

(1) 

where the first (second) term is the Rashba (k-linear 
Dresselhaus) SOI and aju^jj^ is the Rashba (Dresselhaus) 
parameter. The former arises from structure inversion 
asymmetry (SIA) when the potential drop is ori- 

ented preferentially along one direction (say, z at the 
interface between the insulator and the ferromagnet), it 
produces an electric field on the form E w —dzV^,. In the 
presence of finite electron velocity, this field acts on the 
electron spin like an effective magnetic field [IHj . The lat- 
ter form emerges from bulk inversion asymmetry (BIA) in 
strained non-centrosymmetric crystals or interfaces [26|, 
such as ZnS, GaAs etc. The spin-dependent energy shift 
associated with the spin-orbit coupling is on the form 
(H: 



Ej'^,(k) = ±aR{Mxky 



Mykx) ± aoiMxkx - Myky) 



(2) 



when assuming adiabatic spin dynamics (the itinerant 
electron spin is aligned on the local magnetization M). 
When the magnetization is aligned on z, no in-plane spin 
component exist and the energy shift vanishes. But when 
the magnetization lies in the plane of the layer, the en- 
ergy shift is no more zero. As a consequence, due to the 



presence of interfacial SOI, the system depicted in Fig. 1 
is not equivalent when the magnetization is in the plane 
or out-of-the plane. This gives rise to TAMR effects (see 
discussion in Ref. 2^, 23l|). 

In the present work, we consider a SMTJ consisting in 
F/I/N, where F is a ferromagnetic layer, I is an insulator 
and N is a normal metal (see Fig. [T]). In order to capture 
the most important features of the mechanism described 
here, we choose a minimal model only considering the 
most relevant materials' parameters. Matos-Abiague et 
al. [i^ showed that the anisotropic tunneling can be sat- 
isfyingly modeled within the free electron approximation, 
using interfacial Rashba SOI 2J|. The authors showed 
that reasonable agreement with experimental data could 
be obtained with Rashba parameters of the order of 1- 
5eV.A^ [2^. The free electron Hamiltonian of the junc- 
tion reads: 

H = V V + Uz + aRa-(kxz)6(z). (3) 

2 TOr 



Here, an is the Rashba parameter(2J| and rriz is the effec- 
tive mass of the electron, equal to toq in the electrodes 
and nieffmo in the barrier. Uz is the potential of the 
junction, given by: 
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Uo<z<d ^Uo + — - -eV, 



where Uq and d are the barrier height and thickness, 
eV is the bias voltage, J is the s — d exchange cou- 
pling, (T is the vector of Pauli spin matrices and M = 
(sin 6 cos 0, sin 9 sin 0, cos 9) is the magnetization direc- 
tion of the ferromagnetic electrode, as depicted in Fig. 

HI 

The present minimal model satisfyingly reproduces the 
TAMR rcsuhs obtained by Matos-Abiague ct al. [13, 111 
for the case of Fe/GaAs/Au and the one obtained by 
Park et al. [2l[ in (Pt/Co)„/A10x/Pt. In this article, we 
chose the parameters to model Fe/MgO/Au SMTJ, with 
fc],g=1.09nm-\ fc^g=0.4nm-\ fcA„==0.8nm-i, ^^icV 
and d = Inm. As shown in Fig. 2, the TAMR displays 
an angular dependence in cos 29 as already shown in Ref. 

SSI. 



As expressed in Eq. ([T]), the interfacial SOI induces an 
angular momentum transfer between the spin momentum 
and the in-plane orbital momentum of itinerant electrons. 
This mechanism is at the basis of our study. To the first 
order, this transfer is antisymmetric in k and cancels out 
(Hfl(— k) = —HR{k)). However, due to the presence of 
the local s — d exchange, at the second order in SOI, the 
transfer is not zero, yielding an additional spin density 
in the plane of the layers [ux^y in Eq. ([T])]. Interestingly 
this spin density can in turn exert a torque on the local 
magnetization. This can be seen by writing down the 
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spin continuity equation extracted from Eq. (|3]): 

dm . 2J 1 ,r- /^r ^^ / .\ 

— = —V ■ o- ® V - — m X M + — a, Hi^]).{A) 
at m h in 

The first term on the right-hand side is the spatial diver- 
gence of the spin current in the absence of SOI, the second 
term is the torque exerted by the itinerant spin density 
m on the local magnetization M and the last term is 
the torque arising from the angular momentum trans- 
fer driven by the interfacial SOI. In regular spin-valves, 
tunnel junctions or domain walls, the last term usually 
vanishes and the torque is directly related to the spatial 
variation of the spin current (an additional spin relax- 
ation rate is usually inserted in diffusive systems) Q. 
But in the present case, the SOI acts like a source/sink 
for transverse itinerant spin density that competes with 
the spin current. Since the torque is by definition per- 
pendicular to the local magnetization (|M| = 1), it has 
the general form: 



T = T||M X (z X M) +Tj_z x M, 



(5) 



where T\ | and Tj_ are the in-plane and perpendicular com- 
ponents of the torque. In the following, we calculate the 
local itinerant spin density m and integrate this torque 
on the semi-infinite electrodes in order to account for the 
actual torque felt by the magnetization. 



RESULTS AND DISCUSSION 




FIG. 3. (Color online) (a) Influence of Rashba parameter an 
on the equilibrium perpendicular torque T± for different ex- 
changes J; (b) Influence of barrier height U on the equihbrium 
perpendicular torque for different Rashba parameters Or in 
eV.A^ 



in Fig. 3(b): the perpendicular torque slowly decreases 
when increasing the barrier height. 

A number of experimental studies have reported that 
the Rashba coefficient is linear in bias voltage [281 : an = 



+ Q^fl^ (non-linearities appear at larger voltages). 
This linear dependence in semi-conducting 2DEG is as- 
sociated with the modification of the carrier density as a 
function of the bias |28|. In the case of interfacial spin- 
orbit coupling, it is attributed to the variation of the or- 
bital filling at the interface [l^l . Since the perpendicular 
anisotropy is quadratic in a^, we expect a linear depen- 
dence on the bias voltage: Tj, oc « {a^ + 2a^„a\ 



Thick Barriers: Voltage-driven Perpendicular 
Torque 



Thin Barriers: Current-driven In-Plane Torque 



Let us first consider the case of thick tunnel barriers, 
where d > 2nm. This is typi cally the experimental con- 
figuration exploited in Ref. [ill [T^ . We find that the 
equilibrium torque (14=0) is dominated by the perpen- 
dicular component T±, whose angular dependence is re- 
ported on Fig. 2. The angular dependence is on the form 
oc sin 20, which reflects the symmetry of the system: the 
SMTJ is physically equivalent upon the transformation 
9 —9, 9 + IT. Due to cylindrical symmetry, the spin 
torque (and TAMR) docs not depend on (j) (not shown). 
This perpendicular torque can be readily identified with 
the PMA, which competes with the demagnetizing field 
and favors the perpendicular direction. Interestingly, the 
perpendicular torque only weakly depends on (i) the bias 
voltage (the conductance is vanishingly small) and (ii) 
the barrier thickness (not shown - calculated deviations 
are less than 0.1%). 

Fig. 3(a) shows the dependence of the perpendicular 
torque at zero bias as a function of the Rashba param- 
eter for different exchange energies J. Interestingly, the 
torque is quadratic as a function of an. Furthermore, the 
dependence as a function of the barrier height is displayed 



Let us now consider a thin barrier, similar to the ones 
used in regular tunneling spin torque studies Q, d=0.6 
nm. When applying a bias voltage across the junction, 
spin polarized electrons fiow through the tunnel barrier 
and their spin is reoriented due to interfacial SOI, which 
results in an out-of-equilibrium spin torque [isj . In the 
case of thin barriers, one expect that most of the voltage 
drop occurs within the barrier, therefore little modifica- 
tion of the Rashba coefficient is expected. 

In this case, we find that besides the equilibrium per- 
pendicular torque described above, a current-driven in- 
plane torque, Tn, appears. This in-plane torque displays 
a linear dependence on the bias voltage, as shown in 
Fig. 4(a). Interestingly, the out-of-cquilibrium perpen- 
dicular torque is about one order of magnitude smaller 
than the in-plane torque [see inset in Fig. 4(b)]. The 
non-equilibrium in-plane torque rapidly vanishes when 
increasing the barrier height and thickness [Fig. 4(b)]. 

This is an important result since the nature of the spin 
torque itself changes when modifying the barrier thick- 
ness: the non-equilibrium in-plane torque only appears 
in the limit of thin and low barriers, where equilibrium 
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comprising a pinned layer (bottom) and a free layer (top) , 
separated by a thin insulating spacer (generally MgO). 
An additional capping layer, which can be either made 
of heavy metal or very thin insulator (thinner than the 
insulating spacer) , can be inserted on the top of the free 
layer to ensure a significant TAMR-ST. These devices are 
depicted in Fig. [S] In this conventional MTJs, the free 
layer is submitted to a conventional spin torque on the 
form: 



-1.0 -0.5 0.0 0.5 1.0 
Vb(V) 



FIG. 4. (Color online) (a) Influence of bias voltage on the 
non-equilibrium in-plane torque Tn for or G[2,5]eV.A'^; (b) 
Influence of barrier height U on the non-equilibrium in-plane 
torque for d e [0.6, 3]nm. Inset: Influence of barrier height U 
on the non-equilibrium perpendicular torque. In (b), an = 



5eV.A2 and K = O.IV. 



Torque 


STT 


SOI-STT 


TAMR-ST 


Origin 


Magnetic texture 


SOI 


SOI 


Effective field 


Yes in MTJs 


Yes 


Yes, perpendicular 


(anti-)damping 


Yes 


No 


Yes 


Switching 


Yes 


Yes 


Yes 


Precessions 


Yes 


No 


Yes 



TABLE I. Comparison between STT, SOI-ST and TAMR- 
Torque. 



perpendicular torque is always present, giving rise to per- 
pendicular magnetic anisotropy. 



Comparison with STT and SOI-ST 

The present spin torque (referred to as TAMR-ST) can 
be readily compared with the conventional STT [l| and 
the SOI-ST [l4|, as summarized in Table 1. On the first 
hand, it originates directly from the interfacial SOI rather 
than from the inhomogeneous magnetic texture. There- 
fore, it does not need any external polarizer. On the other 
hand, it possesses two components (like STT in MTJs), 
Til and T±, whereas the SOI-ST produces only an effec- 
tive magnetic field [13, [l3|- Furthermore, the in-plane 
component acts like a (anti-)damping (oc M x (z x M)) 
while the perpendicular component (cx: z x M) competes 
with the demagnetizing field. Finally, the angular de- 
pendence of the torque in sin 20 (versus sinO in MTJs) 
is expected to produce original current-driven magneti- 
zation dynamics [18 1. 



DEVICE IMPLICATIONS 

In this section, we discuss the possible implications of 
such a torque in the magnetization dynamics of three 
different devices based on conventional MTJ structures. 



L mtj 



T|7'^M x (P X M) + T""^P X M (6) 



where P is the direction of the polarizer. The presence 
of interfacial SOI at both interfaces of the free layer gen- 
erates an additional TAMR-ST: 

Tsoi = Tp cos 6IM X (z X M) + Ti°' cos 6lz x M (7) 

In Eq. ([7]), we explicitly accounted for the peculiar an- 
gular dependence of the torque in sin 20. Since the per- 
pendicular torque is essentially an equilibrium torque, as 
discussed above, its major influence is to create PMA. 
Most interestingly, the non-equilibrium in-plane torque 
T|^°' is linear in bias voltage, similarly to T™*"', and can 
be in principle of the same order of magnitude. 




ing on the magnetic configure: parallel (T™*-* 



tiparallel (T™*^ 



FIG. 5. Experimental set-ups: (a) perpendicular magnetic 
tunnel junction, (b) in-plane magnetic tunnel junction and 
(c) vortex core. 

The simplest configuration is given in Fig. Ufa), where 
the magnetization of both pinned and free layers are per- 
pendicular to the plane (P = z). Then, both in-plane 
torque will compete against or add to each other, depend- 

0, an- 

r||°') or perpendicular (T™*"') magnetic 
configuration. The dynamics is expected to be much 
more complex in in-plane MTJs [Fig. 5(b)], since P 7^ z. 
Although the TAMR-ST will not participate to the onset 
of the magnetization switching (cos0 « for 9 « 7r/2), it 
will have a strong influence on the magnetization dynam- 
ics itself, favoring either out-of-plane or in-plane states 
depending on the sign of the bias voltage. Finally, Fig. 
Efc ) shows a configuration similar to the one used in Ref. 
[29| , where the spin torque acts on a vortex wall. In the 
absence of TAMR-ST^ the in-plane STT excites strong 
vortex oscillations [29[. In presence of interfacial SOI, 
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the in-plane TAMR-ST is expected to dramatically mod- 
ify the magnetic texture: Since this torque tends to align 
the magnetization either in the plane or out of the plane 
of the layers, it is expected to either broaden or narrow 
the vortex core, depending on the bias voltage direction. 
More detailed will be given elsewhere. 



CONCLUSION 

Spin transport in a SMTJ in the presence of interfacial 
Rashba SOI was investigated theoretically. Besides the 
well-known TAMR effect, we showed that a spin torque 
arises from the presence of interfacial SOI. In the large 
and thick barrier limit, the conductance vanishes and the 
equilibrium torque induces PMA that can be in principle 
controlled by interfacial SOI. In the opposite limit of low 
and thin barriers (large conductance) , a non-equilibrium 
torque arises that allows for magnetization excitations. 
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